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1&  ABSTRACT  |V4x*?N#niO0Marak/ 


An  automated  procedure  has  been  developed  to  perform  seismic  analysis  using  available  historic  and  geologic 
data.  The  objective  of  the  seismicity  study  is  to  determine  the  probability  of  occurrence  of  ground  motion  at  the  site. 
Response  spectra  and  time  history  techniques  are  presented. 
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ZHTROOnCTZOM  >  What  is  a  Seismicity  Study 


The  objective  of  a  seismicity  study  is  to  quantify  the  level 
and  characteristics  of  the  earthquake  ground  motion  which  pose  a 
risk  to  a  site  of  interest.  The  seismicity  study  will  produce  a 
probability  distribution  of  expected  site  acceleration  for  a 
given  exposure  period  and  also  give  an  indication  of  the  frequen> 
cy  content  of  that  motion.  The  approach  taken  in  this  work  is  to 
use  the  historical  epicenter  data  base  in  conjunction  with  geo¬ 
logic  data  where  available  to  best  estimate  the  earthquake  recur¬ 
rence  of  a  region  or  fault.  This  recurrence  relationship  is  used 
to  determine  the  regional  or  fault  contributions  to  the  overall 
site  acceleration  level.  This  becomes  the  basis  for  definition 
of  response  spectra  suitable  for  use  in  structural  design  and 
analysis.  The  procedure  utilizes  three  parts  to  accomplish  the 
study.  The  first  part  creates  a  subset  of  earthquakes  from  the 
general  data  base  and  plots  all  events  within  a  specified  region. 
The  second  part  utilizes  the  epicenter  data  to  perform  a  regional 
analysis  determining  the  magnitude  recurrence  relationship  for 
the  region  which  may  be  adjusted  for  geologic  data  where  known. 
Additionally  the  program  computes  the  probability  of  acceleration 
at  the  site  location  and  gives  plots  of  recurrence  and  accelera¬ 
tion  data.  The  third  part  analyzes  individual  faults.  It  deter¬ 
mines  fault  magnitude  recurrence  and  probability  of  acceleration 
at  the  site  from  an  event  on  each  fault  specified.  Geologic  data 
may  be  used  to  augment  historical  epicenter  date.  Each  part  will 
be  discussed  in  detail  below. 
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System  Requirements 

The  Program  is  designed  to  run  on  standard  desk  top  personal 
computers  using  the  MS  DOS  operating  system  version  3  or  higher. 
The  following  are  required: 

MS  DOS  version  3  or  higher 
640  k  system  memory 
Hard  disk 

80287  math  coprocessor 


optional  devices 


plotter 

printer 


DTIC  QTTALTrf  rFS?r,(Ti’ED  3 


The  following  plotters  are  supported: 


EPSON  FX80  printer  or  compatible  (used  as  plotter) 
Hewlett  Packard  Laser jet  printer  or  compatible 
Hewlett  Packard  Plotters 
Houston  Instruments  Plotters 
Tecktronix  4025 
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Loading  The  Program 

Place  the  first  PROGRAM  DISK  in  Drive  A:  switch  to  Drive  A 
then  type; 


A>  INSTALL  <cr> 

This  will  create  a  directory  on  the  C  hard  disk  named  SEISMIC  and 
load  all  of  the  programs  into  that  directory.  When  complete, 
place  the  next  disk  in  Drive  A  and  type  INSTALL  as  edsove  to 
continue  loading  the  program. 

Place  the  first  DATA  DISK  in  Drive  A:  type  INSTALL  as  above. 
This  will  create  a  directory  named  EPICENTR  and  load  the  data 
files.  Place  the  remaining  DATA  DISKS  in  sequence  in  Drive  A  and 
type  INSTALL  as  above  to  load  all  data. 

Place  the  SAMPLE  DATA  DISK  in  Drive  A:  type  INSTALL  as  edsove. 
This  will  create  a  directory  named  OUTPUT  and  load  the  data 
files. 

See  the  ADVANCED  USER  section  to  deviate  from  these 
locations  for  progrzun  and  data. 

Configuring  The  Program 

The  program  supports  a  number  of  graphic  output  devices;  the 
user  must  specify  which  device  is  being  used  and  how  it  is 
connected.  He  must  also  tell  the  program  where  the  data  files 
are  being  stored.  To  begin  the  configuration  switch  to  the 
program  directory  and  start  the  program  by  typing  the  following: 

C: 

CD\SEISMIC 

SEISMIC 

The  following  screen  will  appear: 


H  OPTIONS 


Directory 

Configure 


USE  AlffiOU  KEYS  THEN  PRESS  D1TER 
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The  program  has  five  choices  which  form  the  opening  menu: 
Options,  Earthquake  Selection,  Regional  Study,  Fault  Study,  and 
Exit.  Use  the  LEFT  and  RIGHT  ARROW  keys  which  are  on  the  number 
pad  keys  with  the  4  and  6  to  move  eunong  the  choices.  Make  sure 
the  HUM  LOCK  key  is  off  to  permit  the  arrow  keys  to  function. 
When  the  OPTIONS  choice  is  selected  a  window  opens  giving  two 
choices:  Directory,  and  CONFIGURE.  Use  the  down  arrow  key  to 
choose  the  CONFIGURE  choice  then  press  ENTER  (or  RETURN) .  The 
following  screen  appears: 


Enter  Plotter  DEVICE  NUnBER 
see  User's  Nanual  e.g.  LPTl  -  1 

pm 

Enter  Plotter  HODEL  NUflBEB 

see  User's  Nanual  e.g.  HP  Lasejet  =68 

Enter  Drive  and  Directory  for  EPICENTER  files 
exanple  CAEQUAKES 

f:veiiiiiakes 

Enter  Drive  and  Directory  for  OUTPUT  files 
exanple  CAOUTPUT 

r:\scisnic 

Any  Gianyes?  V  >'  N 


The  DEVICE  NUMBER  refers  to  the  port  to  which  the  hard  copy 
plotting  device  is  connected;  see  Tedble  1  for  configuration 
options  and  check  the  manual  for  the  hard  copy  plotting  device. 
Enter  the  DEVICE  NUMBER  then  press  ENTER  or  DOWN  ARROW. 

The  MODEL  NUMBER  can  be  obtained  from  Table  2a  for  the 
devices  supported  by  the  program.  Laser  printers  produce  high 
quality  plots  rapidly  and  are  recommended.  Table  2b  gives  a 
matrix  MODEL  NUMBERS  for  various  compatible  printers  which  can 
be  used  to  obtain  plots.  Table  3  gives  the  recommended 
configuration  for  specific  devices.  Enter  the  MODEL  NUMBER  then 
press  ENTER  or  DOWN  ARROW. 

Type  the  directory  name  where  the  epicenter  data  files  are 
located.  If  you  used  the  INSTALL  routine  to  accept  the  default 
directozry  creation  then  type: 

C:\EPICENTR 

Note  the  back-slash  (\)  key  and  the  spelling  of  EPICENTR,  8 
letters  without  the  E.  If  you  chose  to  locate  the  programs  else 
where  enter  the  following: 

Drive  letter: \Directory 
example  D:\EPIC 
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Dttvlc* 


Output  Duvlcu 


Parullul  Port 


PSN:  (PUN:  is  uquivalsnt  to  LFTl:) 
LRl: 

LPn: 

LPT3: 


CON:  Coosolu 


4S00 

4801 

4802 

9600 

9601 

9602 
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T«bl«  2«.  (eeiitlBMd) 


80 

HP  75808,  HP  75858,  or  HP  75868  Drsftlnt 

Plotter  vaiflt  slso  A/A4  to  D/Al  pspor. 

HP  7550A  Grs^ies  Plotter  uslat  else 

A/A4  to  8/A3  pepor. 

HP  7440A  ColorPro  plotter  us inf  siee 

08/A4  pepor. 

85 

HP  75858  or  HP  7586B  Dreftiag  Plotter 
tteint  else  E/AO  p^r. 

90 

Tektronix  4025. 

99 

IBH  color  scabies  oonitor  (CUT). 

Tabl«  2b.  Dot  Itatrix  Friator  Uoogo  by  Ffodol 


*  ■  Tho  prlotor  cn  ooo  this  oodsl  aaobor. 
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Tthim  3.  SttcosMadad  Coaflguratlon 


Ootpat  d«vlc« 

Dsvles 

Modsl 

Epson  FX-80 

0 

5 

Epson  HX*80 

0 

1 

ZBH  Prlntsr 

0 

1 

Csntronics  6LP 

0 

1 

Okldsts  92 

0 

1 

Epson  RX'SO 

0 

1 

Epson  FX*100 

0 

15 

Epson  MX- 100 

0 

11 

LQ«1500 

0 

41 

Okldsts  93 

0 

11 

HI  OHP-Sl 

9600/9650 

51 

HI  DMP-S2 

9600/9650 

51 

HP  7440A 

9600/9650 

80 

HP  7470A 

9600/9650 

20 

HP  747SA 

9600/9650 

30 

HP  7SS0A 

9600/9650 

80 

HP  7S80B 

9600/9650 

80 

HP  7S8SB 

9600/9650 

80/85 

HP  7S86B 

9600/9650 

80/85 

HP  26864 

9600/9650 

60/61 

TsktronlJC  4025 

4800/4850 

90 

ZBH  color  gr^hles 
Monitor 

.99 

99 

No  spaces  are  allowed  in  the  neune.  For  further  information  see 
your  DOS  manual  on  sub-directory  creation  and  naming.  Press  ENTER 
or  DOWN  ARROW  to  advance  to  the  next  item. 

Type  the  directory  name  where  the  output  files  are  located. 
If  you  INSTALLed  the  SAMPLE  DATA  you  created  a  directory  called: 

C:\OUTPUT 

If  you  did  not  then  a  directory  has  not  been  created.  Select  and 
type  a  name  for  the  output  file  location  such  as: 

C:\OUTPUT 

Remember  this  name;  we  will  use  it  again.  Press  ENTER 

Upon  completion,  you  are  asked  if  you  wish  to  make  changes. 
If  everything  is  correct  type  N  for  no  or  else  Y  for  yes  and 
repeat  the  data  entry.  Use  DOWN  ARROW  keys  to  skip  over  accept¬ 
able  answers  and  change  only  what  is  in  error. 

If  you  did  not  use  the  INSTALL  for  SAMPLE  DATA  do  the  following. 
When  N  for  NO  CHANGES  is  pressed  the  program  returns  to  the 
OPTIONS  WINDOW.  Use  the  RIGHT  ARROW  key  to  EXIT.  When  out  of 
the  progreun  create  the  OUTPUT  LOCATION  directory  if  one  does  not 
exist.  Type  the  following: 


CD\ 

MD\OUTPUT 

CD\SEISMIC 

SEISMIC 

The  opening  screen  will  reappear  as  shown  above.  This  time  press 
ENTER  to  accept  the  DIRECTORY  choice.  A  screen  like  the  following 
will  appear  showing  the  earthquake  epicenter  files  on  disk. 


DIRECTORY  F:NEP1CN».EPC 


11  file(s)  found 


Use  cursor  keys  to  select  file  then  press  <ENTER> 
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Press  ENTER  to  continue  and  a  directory  of  the  files  in  the 
output  location  directory  will  appear: 

DiBEcrosy  f;\sdata\».» 


16  file(s)  found 


ONE.PLT 

no. PL 

TCST.3PT 

TUO.PLT 

no.ou 

RESPONSE. PLT 

THREE.PLT 

IEST.30T 

ONE.OUT 

TEST.IOT 

TUO.OUT 

lEST.lPI 

THREE.OUT 

TEST. EOS 

ONE.PL 

Filenane  :  ONE.PLT 


Use  cursor  keys  to  select  file  then  press  <ENIE]i> 


Press  ENTER  to  return  to  the  OPTIONS  window  menu.  We  are  now 
ready  to  begin  a  problem. 

EARTHQOAXB  8ELBCTZON 

The  epicenter  data  base  has  been  prepared  from  the  National 
Oceanic  and  Atmospheric  Administration's  data  base.  Each  file 
covers  a  specified  region  and  contains  date,  latitude,  longitude 
and  magnitude  data.  This  section  is  used  to  create  a  subset  of 
earthquakes  for  a  specific  region  from  the  main  set  of  earthquake 
files.  The  program  searches  a  rectangle  specified  by  maximum  and 
minimum  latitudes  and  longitudes  to  find  all  events  within  the 
box  above  the  specified  magnitude.  To  become  familiar  with  the 
program,  the  user  may  advance  to  the  EQ  SELECTION.  The  follow¬ 
ing  screen  appears: 


9 


EQ  SELECTION 


i.jl]!  a  (ijl  O'  j  I;  I  tj'tKiii  a  1,1 1 V4,  i.,  jl  ;'4liilif.~t!iK‘iiniii 


Specify  Area 
Devise  Data 
Begin  Analysis 


View  Results 
Print  Results 
Plot  Results 
Save  Results 


USE  ARROW  RETS  THEN  PRESS  ENTER 


Use  the  DOWN  ARROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  EQ  SELECTION  window  to  see  the 
output  file  on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS 
will  print  or  plot  the  sample  problem. 

Specify  Data 

This  choice  peirmits  data  entry  for  a  new  problem.  Selecting 
this  choice  will  overwrite  previous  data  examples  unlei.  3  they 
were  saved  to  a  named  file.  You  will  overwrite  the  sample  data 
but  it  can  be  copied  from  the  original  disk  again  if  needed.  The 
following  screen  will  appear: 


Enter  Neximw  Longitiude  (Degrees) 
exenple  120.0 

m 

Enter  Nininun  Longitude  (Degrees) 
example  112.0 

113 

Enter  Naximin  Latitude  (Degrees) 
example  48.0 

4lP;i 

Enter  Hinimum  Latitude  (Degrees) 
example  34.8 

37 

Enter  Site  Longitude  (Degrees) 
example  115.0 

115 

Enter  Site  Latitude  (Degrees) 
example  38.8 

38 

Enter  Hinimum  Nagnitude  Cutoff 
example  3.0 

3.0 

Any  Quinges?  Y  /  N 
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Enter  the  longitude  and  latitude  for  a  rectangle  bounding  the 
problem  and  the  site  longitude  and  latitude.  The  study  area  must 
be  large  enough  to  include  distant  events  which  can  influence  the 
site.  Geologic  data  should  be  consulted  to  look  for  boundaries 
of  tectonic  provinces.  The  site  should  be  near  the  center  of  the 
study  region  unless  otherwise  required  by  a  tectonic  boundary. 

The  region  will  form  the  boundary  for  selecting  events  to  be  used 
in  the  study  and  are  relevant  to  establishing  the  site  seismic 
potential.  For  regional  studies  consideration  should  be  given  to 
selecting  an  area  large  enough  to  establish  the  regional  tectonic 
setting.  Enter  the  MINIMUM  MAGNITUDE  CUTOFF  value,  typically  3.0. 
Events  below  this  level  may  not  have  been  recorded  and  their 
absence  distorts  the  relationship  for  recurrence. 

The  program  computes  the  acceleration  at  the  site  location 
from  each  epicenter  location  where  a  magnitude  is  specified  in 
the  subset  of  earthquakes.  This  computation  uses  an  acceleration 
distance  attenuation  equation.  The  following  attenuation  rela¬ 
tionships  for  acceleration  are  included: 

1.  McGuire  (1978) 

2.  Trifunac  and  Brady  (1975) 

3.  Campbell  (West)  (1982) 

4.  Campbell  (East)  (1982) 

5.  Donovan  and  Bornstein  (1978) 

6.  Joyner  and  Boore  (1981) 

The  user  is  has  the  option  of  selecting  which  to  use.  Since  each 
equation  analyzes  the  data  differently  and  since  there  is  a  large 
uncertainty  in  ground  motion  results  may  differ  depending  upon 
the  equation  selected.  The  user  is  encouraged  to  read  the  refer¬ 
ences.  The  following  screen  is  used  to  enter  the  acceleration 
equation  choice  and  the  hypocenter  depth  if  needed,  usually  10 
miles. 


1  He  Cttire 

2  Trifunac  and  Brady 

3  Canpbell  (Uestem  US) 

4  Canpbell  (Eastern  US) 

5  Donovan  and  Bornstein 

6  Joyner  and  Boore 
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Enter  depth  to  hypocenter  (niles) 


8 


ENTER  CHOICE  FOR  ACCELERATION  EQUATION 

11 


Once  the  acceleration  equation  data  entries  are  completed,  the 
user  is  asked  to  select  epicenter  files  to  be  searched  .  Files 
have  been  divided  into  regions  to  keep  the  search  time  to  a 
minimum.  A  screen  similar  to  the  following  is  shown  from  which 
the  user  may  move  through  the  list  using  the  ARROW  KEYS, 
PAGEUP/DOWN  KEYS  and  then  press  ENTER  to  select  his  choices. 
There  is  no  limit  to  the  number  of  choices. 


Press  ESCAPE  when  done.  The  program  then  shows  the  choices  and 
asks  for  confirmation. 


YOU  HAVE  SELECTED  THE  FOLLOUING  RECORDS: 


7  Calif  1M-1Z3  37-38 

8  Calif  114-124  38-39 

9  Calif  114-125  39-48 

10  Calif  114-125  40-41 

11  Calif  114-125  41-42 


IS  THIS  CORRECT?  Y/N 
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Revise  Data 


If  the  user  wishes  to  revise  a  number  he  may  choose  the 
REVISE  DATA  choice  and  will  be  given  the  same  questions  as  in  the 
SPECIFY  DATA  section  with  the  previous  choices.  Press  enter  or 
TOWN  ARROW  to  accept  the  value  unchanged.  Overwrite  the  revised 
value  completely  to  alter  a  number. 

Begin  Analysis 


This  choice  begins  the  actual  data  search. 
View  Results 


The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 

Print  Results 


PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  a  list  of  epicenters  with  the  computed  site  acceleration  for 
that  event,  and  a  histogram  of  the  distribution  of  acceleration. 

THE  FOLLOWING  IS  AN  EXAMPLE  OUTPUT 


MAX  LONGITUDE  73.000 
MIN  LONGITUDE  69.000 
MAX  LATITUDE  44 . 000 
MIN  LATITUDE  41.000 
MIN  MAGNITUDE  . 000 
MAX  MAGNITUDE  9.000 
SITE  LONGITUDE  72.000 
SITE  LATITUDE  42.000 
ACCELERATION  EQUATION  . 000 


LIST  OF  SITE  EPICENTERS 


YEAR 

LATITUDE 

LONGITUDE 

MB 

MS 

MO 

ML 

AVM 

DISTANCE 

ACCEL 

1976. 

41.66 

69.97 

.00 

.00 

.00 

3.00 

3.00 

106.776 

.003 

1979. 

43.98 

69.80  3 

.80 

.00  4.00 

4.10 

3.97 

177.364 

.004 

1977. 

41.84 

70.70 

.00 

.00  3.10 

.00 

3.10 

67.614 

.006 

1974. 

41.70 

71.50 

.00 

.00 

.00 

2.50 

2.50 

32.976 

.008 

1976. 

41.56 

71.21 

.00 

.00 

.00 

3.50 

3.50 

50.660 

.011 

where 


MB 

MS 

MO 

ML 

AVM 

DISTANCE 

ACCEL 


Body  wave  magnitude 
Surface  wave  magnitude 
Other  magnitude 
Local  magnitude 
Average  magnitude 
Distance  event  to  site 

Site  acceleration  computed  estimate,  g's 
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Plot  File 


This  choice  creates  an  epicenter  plot  for  the  region. 

Figure  1  is  an  exeunple  plot. 

Save  Results 

This  choice  permits  the  user  to  write  the  input,  output  and  plot 
files  to  a  named  file  of  the  user's  choice.  This  prevents  the 
files  from  being  overwritten. 
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73.00 


69.00 


44.00 


41.00 


Note:  numbers  are  approximate  magnitudes 
Figure  1.  Epicenter  plot  for  region. 


15 


REGIONAL  STUDY 


Moving  the  RIGHT  ARROW  key  to  REGIONAL  STUDY  reveals  the 
following  screen: 


' i; J : 'a lill, i! ' ,!  ‘ n )i a-li :) (i' i 1 1 1 'g-  ■"-■■'■'“hI 


Specify  Begion 
Revise  Seta 
Begin  Analysis 


Vieu  Results 
Print  Results 
Plot  Itesults 
Save  Results 
Response  Plot 


USE  ARROU  KEYS  IHQI  PRESS  ENTER 


This  section  is  used  to  conduct  a  regional  seismicity  study 
typically  for  eastern  sites  where  faulting  is  not  known  It  uses 
the  epicenter  subset  created  previously  .  This  section  will 
estimate  the  level  and  characteristics  of  the  earthquake  ground 
motion  which  pose  a  risk  to  the  site  of  interest.  We  will  use 
the  historical  epicenter  data  base  in  conjunction  with  geologic 
data  where  available  to  best  estimate  the  probability  of  site 
acceleration  levels.  This  becomes  the  basis  for  definition  of 
response  spectra  suitable  for  use  in  structural  design  and  analy¬ 
sis. 


For  a  regional  analysis  the  epicenter  data  base  is  used  to 
define  the  regional  magnitude  recurrence  based  on  a  log-linear 
fit  of  the  data  and  the  Richter  A  and  B  recurrence  coefficients 
are  determined. 


Log  (N)  *  A  +  B  M 

The  program  allows  the  user  to  specify  minimum  cutoff  magnitude 
to  enhance  the  fit.  Generally  the  epicenter  data  base  is 
deficient  on  small  events  since  events  less  than  magnitude  3  may 
not  be  large  enough  to  be  recorded  at  distant  seimograph 
stations.  Thus  a  cutoff  of  3  is  usually  used  to  insure  the  fit 
of  the  data  is  through  the  linear  portion  of  the  data.  Figure  2 
illustrates  the  case  in  which  the  initial  data  did  not  use  a 
cutoff  minimum  magnitude.  The  line  should  be  fit  through  the 
linear  portion.  The  program  may  be  repeated,  with  the  user 
specifying  the  A  and  B  coefficients  based  on  his  modification  to 
the  computed  fit  of  the  data. 
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PBOBABILITr  N0T  EXCEEDING  A 


Having  established  the  magnitude  recurrence  relationship,  a 
"floating"  earthquake  analysis  is  performed  using  a  Monte  Carlo 
simulation  routine  A  series  of  events  representing  a  5000  year 
exposure  are  randomly,  spatially  assigned  with  magnitude  shaped 
by  the  recurrence  relationship.  A  maximum  cutoff  magnitude  may 
be  included  to  fit  the  tectonics  of  the  region.  A  probability 
distribution  for  a  specified  exposure  period  of  acceleration  at 
the  site  location  is  computed.  The  simulation  process  is  random 
in  location  and  includes  the  error  associated  with  the  earthquake 
attenuation  equations.  A  histogram  of  magnitude  and  acceleration 
is  computed. 

Use  the  DOWN  ASROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  window  to  see  the  output  file 
on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS  will  print  or 
plot  the  sample  problem. 

Specify  Region 

From  the  REGIONAL  STUDY  window  select  SPECIFY  REGION  and 
press  ENTER.  The  following  questions  will  appear: 


Enter  lUucinuii  Longitinie  (Degrees)  121 

exanple  120.8 

Enter  HiniMin  Longitude  (Degrees)  116 

exanple  112.0 

Enter  Oaxinun  Latitude  (Degrees)  fl 

exanple  dO.O 

Enter  flininun  Latitude  (Degrees)  36 

exanple  3d. 8 

flininun  Earthguake  Hagnitude  Cutoff  3 

exanple  3.0 

flaxinun  Earthguake  flagnitnde  Cutoff  B,S 

exanple  8.0 

Enter  Site  Longitude  (Degrees)  117 

exanple  115.0 

Enter  Site  Latitude  (Degrees)  38 

exanple  38.0 


Ang  Cfianges?  V  ^  N 
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Enter  the  longitude  and  latitude  of  a  rectangle  to  bound  the 
study  area.  This  may  be  a  smaller  than  the  region  chosen  for  the 
EARTHQUAKE  SELECTION  discussed  previously.  Consider  the 
tectonics  of  the  region  in  selecting  bounds  for  the  study.  Be 
sure  to  include  sufficient  distance  from  the  site  so  that  all 
events  which  can  cause  significant  groiind  motion  at  the  site  are 
included.  Enter  the  MINIMUM  MAGNITUDE  CUTOFF,  usually  3  since 
events  less  than  3.0  may  not  have  been  recorded  and  distortion  of 
the  recurrence  estimation  might  result.  Enter  the  MAXIMUM 
MAGNITUDE  of  the  region. 

As  discussed  above,  enter  the  acceleration  equation  and 
hypocenter  depth  if  needed,  usually  10  miles.  Enter  the  EXPOSURE 
PERIOD  of  the  study  in  years. 


1  flc  Cuire 

2  Trifunac  and  Brady 

3  Canpbell  (Uestem  US) 

4  Campbell  (Eastern  US) 

5  Donotwn  and  Bornstein 

6  Joyner  and  Boore 

t 

ENTER  CHOICE  FOR  ACCELERATION  EQUATION 

Enter  depth  to  hypocenter  (niles)  8 

Enter  Exposure  Perion  (Tears)  58 

exanple  50.8 
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At  this  point  you  will  be  asked  whether  to  use  the 
earthquake  data  base  or  to  use  regional  recurrence,  A  and  B 
values.  Enter  Y  to  use  the  epicenter  data.  If  N  is  entered 
enter  the  values  of  A  and  B  as  shown  here: 


Do  you  uish  to  conpute  recurrence. data 
fron  the  epicenter  data  base  (V)es  or  (N)o 
If  you  answer  flo,  you  nust  enter 
Richter  A  and  B  values  for  the  region 
Your  Choice?  Y  ^  N 

Enter  Richter  A  value  li.9 

exanpie  4.2 

Enter  Richter  B  value  -.3 

exanpie  ~.8S 


Revise  Data 

The  user  may  revise  data  once  entered  by  selecting  REVISE 
DATA  from  the  REGIONAL  STUDY  window.  The  same  questions  given  in 
SPECIFY  DATA  are  asked  with  the  previous  responses.  Advance 
through  the  data  using  the  DOWN  ARROW.  Oveirwrite  the  revised 
value  completely. 

Begin  Analysis 

This  choice  begins  the  actual  data  search. 

View  Results 

The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 
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Print  Results 


PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  a  recurrence  data,  the  regional  recurrence  coefficients,  and 
the  site  probability  distribution. 


SAMPLE  OUTPUT 


H/DEG  LONGITUDE  50.906 
M/DEG  LATITUDE  69.057 
MAX  LONGITUDE  73.000 
MIN  LONGITUDE  69.000 
MAX  LATITUDE  44.000 
MIN  LATITUDE  41.000 
MIN  MAGNITUDE  3.000 
SITE  LONGITUDE  71.000 
SITE  LATITUDE  42 . 500 
EXPOSURE  TIME  50.000 
ACCELERATION  EQUATION  . 000 


MAGNITUDE  NUMBER  EVENTS 


.10 

10 

.20 

10 

.30 

10 

.40 

10 

8.70 

0 

8.80 

0 

8.90 

0 

9.00 

0 

YEARS 

COVERED 

RICHTER  FIT  EQ  LOGIO 


GE  M  NUMBER  EVENTS/ YR  GE  M 
.2564 
.2564 
.2564 
.2564 

.0000 

.0000 

.0000 

.0000 

39.00 
)*A  +  B*M 


A  =  3.698  B  =  -.909 


REGIONAL  PROBABILITY  DISTRIBUTION 


CONFIDENCE  ACCELERATION 

MAGNITUDE 

DISTANCE 

LOWER  CONF 

UPPER  CONF 

.9900 

.2276 

5.5186 

13.8088 

.9705 

.9965 

.9802 

.2086 

4.8364 

14.0230 

.9537 

.9918 

.9704 

.1973 

5.2427 

16.7578 

.9390 

.9864 

.9608 

.1796 

5.6075 

20.4107 

.9254 

.9805 

.9512 

.1500 

4.9374 

13.8097 

.9125 

.9742 

.9418 

.1461 

5.0632 

13.8910 

.9002 

.9677 

.9324 

.1425 

5.4726 

13.8134 

.8883 

.9610 

.9231 

.1397 

6.0602 

14.2851 

.8768 

.9541 

.9139 

.1371 

5.0184 

13.8089 

.8656 

.9485 

.9048 

.1356 

5.2686 

13.9865 

.8547 

.9402 

.8958 

.1301 

4.2273 

14.7429 

.8440 

.9331 

.8869 

.1300 

4.7795 

16.7417 

.8335 

.9260 

.8781 

.1282 

4.6073 

14.9659 

.8233 

.9188 
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.1262  5.2475  15.1447  .8133  .9117 
.1239  5.5734  15.8349  .8034  .9044 
.1237  5.3539  18.6373  .7939  .8972 

more 


HISTOGRAMS  COVERING  5000.0  YEARS 


NO.  P(X) 

F(X) 

HISTOGRAM  MAGNITUDE 

86  .086 

.086 

3.307 

lilthlilUtliiUtlHilUtlfk*  It  *****11**1111  It  it  It  **11  It  *********** 

109  .109 

.195 

3.450 

********************************************* 

73  .073 

.268 

3.593 

********************************************* 

77  .077 

.345 

3.736 

************************************' 

more 

************************************ 

68  .068 

.413 

3.879 

51  .051 

.464 

4.022 

********************************************* 

66  .066 

.530 

4.165 

********************************************* 

33  .033 

.563 

4.308 

********************************* 

34  .034 

.597 

4.451 

********************************** 

36  .036 

.633 

4.594 

************************************ 

35  .035 

.668 

4.737 

*********************************** 

36  .036 

.704 

4.880 

************************************ 

33  .033 

.737 

5.023 

********************************* 

31  .031 

.768 

5.166 

******************************* 

20  .020 

.788 

5.309 

******************** 

NO.  P(X) 

F(X) 

more 

HISTOGRAM  ACCELERATION 

563  .563 

.563 

.007 

********************************************* 

214  .214 

.777 

.021 

************************************  more 

86  .086 

.863 

.036 

************************************* 

47  .047 

.910 

.050 

********************************************* 

18  .018 

.928 

.064 

****************** 

13  .013 

.941 

.079 

************* 

9  .009 

.950 

.093 

********* 

9  .009 

.959 

.107 

********* 

6  .006 

.965 

.121 

****** 

6  .006 

.971 

.136 

****** 

4  .004 

.975 

.150 

**** 

2  .002 

.977 

.164 

** 

3  .003 

.980 

.179 

*** 

2  .002 

.982 

.193 

** 

2  .002 

.984 

.207 

** 

1  .001 

.985 

.221 

* 

0  .000 

2  .002 

.985 

.987 

.236 

.250 

** 

more 


Plot  File 

This  choice  creates  a  plot  of  recurrence  and  site 
acceleration  probability  for  the  region.  Figure  2  is  an  example 
plot. 


.8694 

.8607 

.8521 
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Save  Results 


This  choice  permits  the  user  to  write  the  input,  output  and 
plot  files  to  a  named  file  of  the  user's  choice.  This  prevents 
the  files  from  being  overwritten. 

Response  Plot 

This  choice  creates  a  standard  shaped  site  independent 
response  plot  on  your  plotting  device.  Figure  3.  The  user  is 
requested  to  provide  the  base  ground  motion  level  in  g's. 
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VELOCITT  (IN/SEC) 


1000.00 


DAMPING  2,  5,  7.  AND  10  PERCENT 


.01  .10  1.00  10.00  100.00 

PERIOD  (SEC) 


Figure  3.  Typical  response  spectra. 
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FAULTS  STUDY 


For  western  sites  where  fault  locations  are  known  in  more 
detail,  a  study  can  be  performed  by  specifying  the  location  of  a 
fault  in  terms  of  coordinates  of  several  points  defining  line 
segments.  All  earthquakes  within  a  specified  distance  or 
boundary  from  the  fault  line  under  study  are  made  a  subset,  and 
the  recurrence  of  the  fault  is  calculated.  A  probability 
analysis  is  then  performed  to  calculate  expected  site 
acceleration  and  causative  earthquake  magnitude  and  epicentral 
distance  for  that  fault.  The  program  randomly  selects  the 
epicenter  of  an  earthquake  somewhere  along  the  specified  length 
of  the  fault.  Using  the  fault  recurrence  data  in  terms  of 
Richter  coefficients  and  maximum  earthquake  magnitude  associated 
with  the  fault,  the  program  determines  and  earthquake  magnitude 
and  the  length  of  fault  break  (assumed  centered  on  the 
epicenter) ,  and  then  calculates  the  distance  of  the  site  to  the 
fault  break  (the  epicentral  distance  and  the  hypocentral 
distance) .  These  distances,  along  with  the  acceleration- 
magnitude  attenuation  relationship  with  its  uncertainty  defined 
by  the  standard  deviation  for  that  level  of  motion,  give  the  site 
acceleration.  The  process  is  repeated  using  a  Monte  Carlo  scheme 
to  produce  a  list  of  site  accelerations  and  related  causative 
magnitude  and  epicentral  distance  thus  defining  the  site's 
probability  distribution.  It  is  important  to  note  that  the 
program  is  random  for  each  and  every  fault  in  the  following: 

a)  Location  along  fault  length,  2  dimensional 

b)  Magnitude  shaped  by  recurrence  coefficient 

c)  Acceleration  level  using  mean  and  standard  deviation 
relationship  of  magnitude  -  distance 

Provisions  are  included  to  use  recurrence  data  from  slip 
analysis  or  regional  seismicity  in  lieu  of  fault  specific  data. 

The  program  determines,  tabulates,  and  plots  the  probability  of 
not  exceeding  various  levels  of  acceleration  in  the  time  period 
specified.  These  data  are  available  for  all  individual  faults,  and 
then  are  combined  for  all  faults  acting  together.  The 
determination  of  total  risk  to  the  site  is  of  importance  for 
establishing  design  levels. 

A  single  recent  large  event  may  release  strain  built  up  over 
hundreds  or  thousands  of  years.  As  such,  it  might  indicate  a 
period  of  less  activity  in  the  immediate  future.  However,  since 
the  data  base  is  relatively  short,  the  return  time  for  this  event 
might  be  erroneously  indicated  as  much  less.  For  example,  if 
this  were  a  500-year  event  and  occurred  during  a  50-year  data 
base  its  return  might  be  estimated  at  0.02  rather  than  0.002. 

The  plotted  data  points  and  line  of  best  fit  are  determined  by 
the  computer  analysis  using  regression  analysis  techniques. 

These  should  be  reviewed  and  judgment  used  to  adjust  this  type  of 
datxim  point  that  will  clearly  plot  significantly  higher  than  the 
linear  portion  of  the  recurrence  data. 

The  procedure  is  intended  to  be  repeated  several  times, 
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during  which  the  engineer  can  compare  geologic  data  and  lines  of 
best  fit  from  historic  data  and  converge  on  the  best  estimate 
using  his  judgment. 

The  FAULTS  STUDY  window  reveals  the  following  choices: 


FAULTS  STUDy 


1 0^9090911 

Revise  Data 
Begin  Analysis 


View  Results 
Print  Results 
Plot  Results 
Plot  Faults 
Save  Results 
Response  Plot 


Edit  Faults. Lst 


USE  ARROU  KEYS  THEN  PRESS  ENTER 


Use  the  DOWN  ARROW  key  and  then  choose  REVISE  DATA  to  see  the 
example  problem.  Press  RETURN  to  accept  each  value  unchanged. 

Then  select  VIEW  RESULTS  from  the  window  to  see  the  output  file 
on  screen.  Selecting  PRINT  RESULTS  or  PLOT  RESULTS  will  print  or 
plot  the  sample  problem. 

Specify  Data 

This  is  the  data  entry  section  for  a  new  problem.  Selecting 
this  choice  shows  the  following  screen: 
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Enter  naxinuR  Lonsitiude  (Degrees) 
exmple  12B.0 


Enter  tlininim  Longitude  (Degrees)  ISliii; 

example  112.0 

Enter  llaxiMM  Latitude  (Degrees)  Hlf 

example  40.0 

Enter  llinimum  Latitude  (Degrees)  piiiP 

example  34.0 

ninimum  ttagnitude  Cutoff  liilii 

example  3.0 

Enter  Site  Longitude  (Degrees)  11?!  1; 

example  115.0 

Enter  Site  Latitude  (Degrees) 
example  38.0 

Enter  Exposure  Period  (Pears)  HlIB 

example  58.0 

Ang  ChamieeT  V  /  N 

Enter  the  longitude  and  latitude  of  a  rectangle  to  bound  the 
study  area.  This  may  be  a  smaller  than  the  region  chosen  for  the 
EARTHQUAKE  SELECTION  discussed  previously.  Consider  the 
tectonics  of  the  region  in  selecting  bounds  for  the  study.  Be 
sure  to  include  sufficient  distance  from  the  site  so  that  all 
events  which  can  cause  significant  ground  motion  at  the  site  are 
included.  Enter  the  MININUM  MAGNITUDE  CUTOFF,  usually  3  since 
events  less  than  3.0  may  not  have  been  recorded  and  distortion  of 
the  recurrence  estimation  might  result.  Enter  the  EXPOSURE 
PERIOD  of  the  study  in  years. 
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Sn'teir  ^he  acceleiration  e(jua^ion  use  as  discussed  above. 

1  He  Caire 

2  Trifttnac  and  Brady 

3  Canpbell  (Uesiem  US) 

4  Canpbeli  (Eastern  US) 

5  Donovan  and  Bomstein 

6  Joyner  and  Boore 

■ 


ENTER  CHOICX  FOR  ACCELEMTION  EQUATION 


The  user  is  then  given  a  menu  of  known  faults  in  the  California 
area.  This  menu  is  NOT  complete;  but  is  meant  as  a  vehicle  to 
reduce  data  entry.  The  user  may  add  to  this  list;  this  will  be 
discussed  later. 


SAN  (UHEHTE 

PALOS  OERDE 

NEUPORT'INGIEUOOD 

ROSE  CANYON 

UNNANED 

ELSINORE 

WITTIER 

SAN  JACINTO 

COYOTE  CREEX 

SUPERSTITION  NT 

SUPERSTITION  HILL 

INPERIAL 

BANNING 

SAN  ANDREAS 

XT  LOCAL 

CRISTIANITOS 

ALISO 


Use  the  ARROW  KEYS  or  press  the  first  letter  of  the  fault  name 
desired  to  move  throught  the  list.  Press  ENTER  to  select  a  fault* 
up  to  IS  faults  may  be  selected. 
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The  question  is  then  asked  whether  to  enter  additional  faults. 
If  the  answer  is  yes  the  following  screen  is  shown: 


Enter  Fault  Nane 
exanple  San  Andreas 

Enter  Haxinun  Hagnitude 
exanple  7.0 


Enter  Characteristic  Magnitude  Hi 

exanple  7.5 


Enter  leturn  Tine  in  Years 
exanple  250.0 


Enter  Hesignator  1  or  2  or  3  IS 

1  for  2  line  segnent 

2  for  hox  segnent 

3  for  specification  of  Richter  A  and  0 


Enter  llepth  to  Hgpocenter  (niles)  81 

exanple  5.0 


Ang  Changes?  Y  /  N 


Enter  the  name  of  the  fault  and  the  faults  MAXIMUM  MAGNITUDE. 

The  CHARACTERISTIC  MAGNITUDE  is  the  magnitude  which  a  geological 
evidence  shows  has  a  frequent  history  of  occur ing  as  a  major 
event  and  the  RETURN  TIME  requested  is  the  return  time  of  that 
event.  This  event  is  added  to  the  seismicity  computed  by  the 
historical  data.  The  CHARACTERISTIC  MAGNITUDE  may  exceed  the 
MAXIMUM  MAGNITUDE  Which  is  the  cutoff  for  the  recurrence 
calculation  based  on  the  historical  data  or  the  input 
coefficients. 
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The  fault  may  be  described  by  three  choices: 

1.  A  2  line  segment  where  events  a  specified  distance  from 
the  fault's  line  segments  are  Included  In  the  subset  of 
earthquake  events  used  to  calculate  the  recurrence  of  the  fault. 

2.  A  2  line  segment  where  earthquake  events  within  a  4  sided 
region  are  used  to  calculate  the  recurrence  of  the  fault. 

3.  A  2  line  segment  where  the  recurrence  of  the  fault  Is 
specified  In  terms  of  the  Richter  A  and  B  values. 

Enter  your  choice  for  the  specific  fault  being  defined. 

If  you  select  1: 

The  following  questions  will  be  asked 


The  following  questions  will  be  asked 


Enter  Paint  1  Longitia4e  (Segrees) 
emnple  12B.0 

120:13 

Enter  Point  1  Latitude  (Degrees) 
exanple  46.8 

<8 

Enter  Point  Z  Longitude  (Degrees) 
exanple  118.8 

118 

Enter  Point  2  Latitude  (Degrees) 
exanple  38.8 

37.5 

Enter  Point  3  Longitude  (Degrees) 
exanple  117.8 

116.3 

Enter  Point  3  Latitude  (Degrees) 
exanple  37.8 

37 

Distance  fron  fault  to  include  Events  (niles) 
exanple  18.8 

10 

Any  Changes?  Y  /  N 
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These  define  the  line  segments  and  distance  away  from  the  fault, 
see  Figure  4. 

If  you  select  2: 


The  following  questions  will  be  asked  as  above. 


Enter  Point  1  Lonjitinie  (Degrees) 
exsnple  128.8 

12811 

Enter  Point  1  latitude  (Degrees) 
exanple  48.8 

48iiP 

Enter  Point  2  Longitude  (Degrees) 
exanple  118.8 

HB 

Enter  Point  2  Latitude  (Degrm) 
exanple  38.8 

MB 

Enter  Point  3  Longitude  (Degrees) 
exanple  117.8 

mm 

Enter  Point  3  Latitude  (Degrees) 
exanple  37.8 

371. 

Additionally  coordinates  for  a  4  sided  region  must  be  entered. 
Start  with  the  uppermost  right  point  and  proceed  CLOCKWISE.  See 
Figure  5. 


EniER  4  POINTS  FOR  lOX  STARTINC  AT  UPPER  RICH!  AND  CO  OOCRUISE 


Enter  Point  1  Longitiude  (Degrees) 
exanple  128.8 

188 

Enter  Point  1  Latitude  (Degrees) 
exanple  46.8 

Enter  Point  2  Longitude  (Degrees) 
exanple  118.8 

118.5 

Enter  Point  2  Latitude  (Degrees) 
exanple  39.5 

38 

Enter  Point  3  Longitude  (Degrees) 
exanple  118.5 

128.5 

Enter  Point  3  Latitude  (Degrees) 
exanple  37.8 

37.5 

Enter  Point  4  Longitude  (Degrees) 
exanple  128.5 

128 

Enter  Point  4  Latitude  (Degrees) 
exanple  37.5 

48 
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If  you  select  3: 


Enter  the  points  for  the  2  line  segments  as  eOsove. 


Enter  Point  1  Longitiudle  (Degrees) 

128 

exanple  128.8 

Enter  Point  1  Latitude  (Degrees) 

48 

exanple  48.0 

Enter  Point  2  Longitude  (Degrees) 

118 

exanple  118.8 

Enter  Point  2  Latitude  (Degrees) 

37;5 

exanple  38.0 

Enter  Point  3  Longitude  (Degrees) 

116.3 

exanple  117.0 

Enter  Point  3  Latitude  (Degrees) 

37 

exanple  37.8 

ftny  Changes?  V  /  N 


Then  enter  the  A  and  B  values. 


Enter  Richter  A  for  Fault 
exanple  4.8 

Enter  Richter  B  for  Fault  '"P  i 

exanple  -.89 
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Revise  Data 


The  user  may  revise  data  once  entered  by  selecting  REVISE 
DATA  from  the  FAULTS  STUDY  window.  The  same  questions  given  in 
SPECIFY  DATA  are  asked  with  the  previous  responses.  The  user  may 
revise  the  data  specified  in  the  predefined  faults  list. 


N  VOU  UISH  TO  REVISE  DATA  FOR  THE  FAULTS 
you  SELECTED  FROH  THE  PREDEFINED  LIST 
(VIES  /  (NIO 


Selection  of  N  for  NO  allows  the  user  to  accept  or  revise 
the  selection  of  faults  from  the  predefined  list  of  faults  in  the 
menu.  The  data  for  these  faults  can  not  be  revised  for  this 
choice.  However,  the  user  may  revise  the  data  for  the  faults  he 
entered.  Selection  of  Y  for  YES  treats  the  predefined  faults  as 
if  they  were  user  entered  and  allows  the  user  complete  freedom  to 
change  all  values.  In  editing  faults  the  user  may  skip  a  fault 
by  pressing  the  ESCAPE  KEY  after  the  fault  name  appears  rather 
than  the  ENTER  or  DOWN  ARROW  KEY.  This  advances  to  the  next 
fault  leaving  the  previous  fault's  values  unchanged. 

NOTE  AND  WARNING 

THE  PREDEFINED  FAULTS  ARE  REASONABLE  FIRST  ESTIMATES.  THEY 
WILL  NOT  SUIT  ALL  CASE  STUDIES.  IN  PARTICULAR,  THE  USER  SHOULD 
PAY  ATTENTION  TO  THE  DEFINITION  OF  THE  REGION  SURROUNDING  THE 
FAULT  LINE  TO  INCLUDE  EARTHQUAKE  EPICENTERS  FOR  THE  FAULT 
RECURRENCE  RELATIONSHIP.  THE  INTENT  IS  TO  INCLUDE  ONLY  THOSE 
EVENTS  WHICH  CAN  BE  ATTRIBUTED  TO  THAT  FAULT  AND  EXCLUDE  THOSE 
FROM  OTHER  FAULTS.  BASED  ON  YOUR  SPECIFIC  PROBLEM  IT  MAY  BE 
NECESSARY  TO  CHANGE  THE  SELECTION  NODE  FROM  A  STANDARD  DISTANCE 
AWAY  FROM  THE  FAULT  TO  A  QUADRILATERAL  DEFINITION.  AFTER  A 
PRELIMINARY  ANALYSIS  FIRST  RUN  SOUTTION,  THE  USER  SHOULD  EXAMINE 
EACH  RECURRENCE  CURVE  AND  INCORPORATE  GEOLOGIC  DATA.  THE  SLOPE  OF 
THE  RECURRENCE  LINE  SHOULD  BE  CHECKED  TO  INSURE  IT  PASSES  THROUGH 
THE  LINEAR  PART  OF  THE  EVENTS  DATA.  THE  USER  SHOULD  REVISE  THE 
FAULT  DATA  FOR  THE  ALPHA  AND  BETA  COEFFICIENTS  COMPUTED  FROM  THE 
INITIAL  PRELIMINARY  ANALYSIS.  THIS  ITERATION  CONTROLS  THE 
QUALITY  OF  THE  ANALYSIS.  THE  ACCURACY  OF  THE  ANALYSIS  IS  A 
FUNCTION  OF  THE  EFFORT  SPENT  BY  THE  USER  TO  DEVELOP  THE  MODEL. 

THE  PROGRAM  HAS  THE  CAPABILITY  OF  PRODUCING  HIGHLY  ACCURATE 
RESULTS. 
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Begin  Analysis 

This  choice  begins  the  actual  data  search. 

View  Results 

The  VIEW  RESULTS  choice  permits  the  user  to  see  the  output 
file  from  an  analysis  on  the  screen.  It  must  be  run  after  an 
analysis  has  been  performed  and  data  exists. 


Print  Results 

PRINT  RESULTS  prints  the  output  files.  The  output  consists 
of  fault  recurrence  data,  the  fault  recurrence  coefficients, 
and  the  fault  and  total  site  probability  distribution. 

SAMPLE  OUTPUT 


M/DEG  LONGITUDE 

57.710 

M/DEG  LATITUDE 

69.057 

MAX  LONGITUDE 

119.000 

MIN  LONGITUDE 

115.000 

MAX  LATITUDE 

34.500 

MIN  LATITUDE 

32.000 

MIN  MAGNITUDE 

3.000 

SITE  LONGITUDE 

117.360 

SITE  LATITUDE 

33.300 

EXPOSURE  TIME 

50.000 

ACCELERATION  EQUATION 
SAN  CLEMENTE 

.000 

INDIVIDUAL  FAULT  STUDY 

FAULT  COORDINATES 

118.700  32.200 

118.300  32.800 

117.800  32.650 


FAULT  MAX  CREDIBLE  EARTHQUAKE  7.70 


FAULT  EPICENTERS 

LONGITUDE 

LATITUDE 

MAGNITUDE 

DIST 

ACC 

117.800 

32.583 

4.500 

55.645 

.018 

117.833 

32.800 

4.000 

44.015 

.017 

117.833 

32.716 

4.400 

48.699 

.020 

117.866 

32.800 

3.100 

more 

45.221 

.009 

AVE  ACCELERATION  .0110 

MAX  ACCELERATION  .0378 


37 


EARTHQUAKE  RECURRENCE 


MAGNITUDE 

NUMBER 

EVENTS 

GE  M 

NUMBER  EVENTS/ YR  GE  M 

.10 

45 

.6164 

.20 

45 

.6164 

.30 

45 

.6164 

.40 

45 

.6164 

8.70 

0 

.0000 

8.80 

0 

.0000 

8.90 

0 

.0000 

9.00 

0 

.0000 

YEARS 

COVERED 

73.00 

RICHTER  EQUATION  LOGIO(N)  =  A  +B  M 
A  =  1.640  -.601 


FAULT  REGION  ACCELERATION  PROBABILITY 


CONFIDENCE  ACCELERATION  MAGNITUDE  DISTANCE  LOWE  R  CONF  UPPER  CONF 


.9900 

.9802 

.9704 

.9608 

.9512 

.9418 

.9324 

.9231 

.9139 

.9048 

.8958 

.8869 


.3014 

.2364 

.2084 

.1984 

.1799 

.1795 

.1606 

.1571 

.1548 

.1306 

.1292 

.1287 


7.2624 

6.7600 

6.0210 

7.3362 

6.7060 

5.9095 

6.6048 

7.4193 

6.9925 

6.7041 

7.5633 

5.6727 


30.0969 

30.0901 

34.3689 

30.2334 

30.1197 

30.5953 

53.9805 

55.4437 

31.2901 

34.9663 

30.1013 

30.0893 


,9705 

,9537 

.9390 

,9254 

.9125 

.9002 

.8883 

.8768 

.8656 

.8547 

.8440 

.8335 


.9965 

.9918 

.9864 

.9805 

.9742 

.9677 

.9610 

.9541 

.9485 

.9402 

.9331 

.9260 


A  SIMILAR  DATA  TABULATION  IS  GIVEN 
FOR  THE  TOTAL  PROBABILITY  ACCELERATION 
FROM  ALL  OF  THE  FAULTS 


Plot  File 

This  choice  creates  a  plot  of  recurrence  and 
probability  data  for  each  fault  and  the  total  site  probability  of 
acceleration.  Figure  6  is  an  example  plot. 
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Geologic  Data 


Original  line 


MRGNITUDE 


Note:  The  user  should  review  the  fit  of  the  recurrence  data 
shown  by  the  line.  This  line  should  pass  through  the  linear 
portion  of  the  data.  Sometines  when  the  niimber  of  events  is 
limited  or  the  minimum  event  is  not  specified  the  line  may 
miss  that  region.  The  user  should  then  rerun  the  problem 
revising  the  fault  data  by  entering  the  A  and  B  values  from 
a  new  fit  correctly  through  the  data. 


Figure  6a.  Recurrence  for  a  fault. 


PMBABILITr  N0T  EXCEEDING  A 


Plot  Faults 


This  choice  creates  a  plot  of  the  faults  entered  in  the 
SPECIFY  DATA  or  REVISE  DATA  option.  The  epicenters  in  the 
EQUAKES  file  are  plotted  also.  This  choice  requires  that  the 
data  be  specified  for  the  faults  and  the  epicenter  file  be 
created.  However  the  analysis  need  not  be  performed,  so  this  may 
be  used  to  check  the  events  near  faults. 

The  plotting  procedure  draws  the  grid  of  the  region 
specified  and  then  draws  the  faults  and  plots  the  epicenters. 
Often  the  fault  extends  beyond  the  region  specified.  IT  IS  NORMAL 
IN  SUCH  CASES  TO  SEE  AN  ERROR  MESSAGE  STATING  UNPLOTTED  VECTORS 
OR  CLIPPED  VECTORS  for  the  elements  which  go  beyond  the  plot 
boundaries . 


Save  Results 

This  choice  permits  the  user  to  write  the  input,  output  and 
plot  files  to  a  named  file  of  the  user's  choice.  This  prevents 
the  files  from  being  overwritten. 

Response  Plot 

This  choice  creates  a  standard  shaped  site  independent 
response  plot  on  your  plotting  device.  Figure  3.  The  user  is 
requested  to  provide  the  base  grotuid  motion  level  in  g's. 

Edit  Faults . Lst 

This  option  permits  the  user  to  change  the  values  of  the 
data  used  to  define  the  faults  in  the  predefined  faults  list 
which  appears  as  a  menu  from  which  the  user  may  select  faults  for 
an  analysis.  Appendix  A  explains  how  to  use  this  option  and  also 
explains  how  to  calculate  the  recurrence  coefficients,  A  and  B 
from  a  recurrence  plot.  The  appendix  gives  a  list  of  faults  and 
the  recurrence  plots  for  all  the  faults  defined  in  the  FAULTS. LST 
data  file  of  predefined  faults.  The  user  may  add  to  this  list. 

NOTE  recurrence  plots  for  faults  with  predefined  recurrence 
coefficients  are  not  plotted  to  save  time;  refer  to  Appendix  A 
for  the  plots. 


42 


Discussion  of  Faults 

Seismometers  have  been  installed  near  known  active  faults  to 
record  microearthguakes .  The  events  recorded  range  in  magnitude 
from  0.5  to  1.5  and  closely  trace  the  faults.  There  are  also 
regions  where  few  microearthguakes  occurred.  The  San  Andreas 
fault  north  of  the  Sargent  fault  exhibits  less  activity  than 
southern  portions,  perhaps  indicating  it  is  locked.  However, 
sufficient  microearthguakes  have  occurred  to  show  the 
continuation  of  the  fault.  It  may  be  reasoned  that  a  small 
failure  might  occur  before  a  major  rupture  occurs ; alternatively , 
a  large  nimber  of  microearthguakes  demonstrate  active  creep  that 
may  be  sufficient  to  prevent  sizedale  strain  accumulation  and 
preclude  a  large  event. 

Since  earthguakes  are  associated  with  faults,  it  might  be 
thought  that  epicenters  should  precisely  overlay  the  fault 
location.  This  is  not  the  case,  because  the  distribution  of 
seismometers  is  uneven  and  has  changed  with  time.  There  are 
limitations  in  the  accuracy  of  the  technigues  used  to  locate 
epicenters,  principally  from  variations  in  assumed  propagation 
velocities.  Further  explanation  for  the  location  of  epicenters 
being  off  their  associated  fault  comes  from  the  simplified  model 
used  to  locate  them.  The  center  of  earthguake  energy  is  located 
at  the  focus.  For  an  inclined  fault,  the  surface  location 
(epicenter)  is  a  distance  removed  from  the  surface  fault 
location.  It  is  only  in  vertical  faults  that  one  might  expect 
the  epicenter  to  lie  on  the  fault. 

Krinitzsky  (1974)  concludes  that  earthguakes  can  be  related 
to  existing  faults  and  that  the  .possibility  of  formation  of  new 
faults  should  not  be  considered  in  design.  Large  earthguakes 
reguire  fault  breaks  of  considerable  distance.  The  uncertainties 
associated  with  earthguakes  generated  from  faults  not  previously 
known  can  occur  only  for  small  events.  Generally  in  the  western 
United  States,  the  extent  of  geologic  investigation  precludes  a 
large  fault  from  remaining  unknown.  However,  there  are 
uncertainties  associated  with  eastern  earthguakes.  For  example, 
causative  faults  responsible  for  the  New  Madrid  earthguakes  of 
1811  and  1812  have  not  yet  been  identified.  This  may  be  the 
result  of  insufficient  geologic  investigations.  The  importance 
of  considering  the  extent  and  guality  of  geologic  investigations 
is  evident . 

A  period  of  demonstrated  guiescence  over  a  geological  time 
period  indicates  inactivity  of  the  fault  and  probable  continued 
inactivity.  However,  inactivity  over  a  period  of  historic 
recording  (50  to  100  years)  does  not  imply  future  inactivity. 
Rather  it  may  point  to  a  "locked"  region  through  which  a  major 
fault  rupture  may  propagate.  Two  earthguakes  to  produce  damage 
in  southern  California  (Arvin-Tehachapi  1952;  San  Fernando,  1971) 
occurred  on  faults  lacking  major  historic  activity.  With  the 
exception  of  the  San  Jacinto  fault  system  every  Icnown  event 
greater  than  magnitude  6  in  southexm  California  has  occurred  on  a 
fault  without  prior  major  historic  activity.  It  must  be 
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recognized  that  the  accuracy  of  an  incomplete  data  base  is  very 
limited  when  extrapolated  for  return  period  greatly  exceeding  the 
length  of  the  period  of  recorded  data.  Furthermore,  aftershocks 
must  be  distinguished  from  main  shocks.  An  area  having  recently 
undergone  a  large  event,  which  releases  strain  built  up  for 
hundreds  or  thousands  of  years,  is  probably  safe  against  that 
size  release  in  the  near  future.  Thus  a  recent  event  on  a  fault 
might  actually  indicate  safety  in  the  immediate  future  rather 
than  an  indication  of  activity.  A  single  event  by  itself  cannot 
give  an  accurate  measure  of  return  time.  The  limitations  of  the 
historic  data  base  can  be  reduced  by  considering  the  geologic 
data  available  to  assist  in  the  definition  of  fault  activity. 

Faults  are  classified  into  several  types  according-  to  the 
relative  displacement  of  the  sides  of  the  fault.  Three  principal 
types  are  normal,  reverse,  and  strike-slip.  The  main  fault  may 
appear  as  a  single  break  or  as  a  parallel  series  of  breaks. 
Surface  rupture  may  or  may  not  occur.  Variations  in  displacement 
occur  along  the  fault.  Some  fault  displacement  occurs 
continuously  at  a  very  slow  rate  and  is  termed  as  "creep".  Creep 
rates  have  been  established  for  some  active  faults  in  California. 
Estimates  of  earthquake  recurrence  can  be  sometimes  made  by 
observation  of  fault  displacements  of  exposed  layers,  in  a  trench. 
From  age  dating  it  is  possible  to  identify  the  age  of  various 
strata  and  also  to  identify  layer  displacements. 

ASSESSMEMT  OF  GROUND  MOTION  ATTENUATION  REXATZON8HZF8 

A  Study  was  perfomned  evaluating  the  significance  and 
difference  of  available  earthquake  attenuation  relationships. 
Data,  although  apparently  plentiful,  are  actually  very  limited. 
Events  are  recorded  by  strong  motion  accelerographs  often  located 
within  structures.  The  effect  of  the  response  of  the  structure 
often  influences  the  reading.  The  assessment  of  the  location  and 
separation  distance  of  many  early  recorded  events  is  of 
questionable  accuracy.  Many  researchers  create  subsets  of  the 
data  often  defining  their  own  meanings  to  distance  and  magnitude 
and  often  limit  the  inclusion  of  records  to  specific  conditions 
such  as  rock  sites  or  events  of  only  a  limited  magnitude  range. 
Campbell  (1982)  for  example  computes  in  his  relationship  the 
average  acceleration  of  the  two  horizontal  directions.  Campbell 
states  that  the  maximum  direction  value  should  be  1.13  times  the 
computed  value.  Joyner  and  Boore  (1981)  use  the  seismic  moment 
as  the  basis  for  magnitude  determination  while  Campbell  uses 
local  magnitude  for  events  less  than  6.0  and  surface  wave 
magnitude  for  events  greater  than  6.0.  Each  investigator  uses 
judgment  and  perception  of  the  elements  incorporated  into  the 
formulation.  This  subjective  approach  makes  seismic  attenuation 
prediction  more  of  an  art  since  the  scatter  in  the  data  is  large. 
Recent  work  suggests  that  magnitude  may  not  have  the  same 
significance  for  ground  motion  levels  close  to  a  fault  as  it  does 
at  greater  distances.  The  analytical  procedures  generally  assume 
motion  to  be  related  to  magnitude  and  distance  as  independent 
variables.  However,  it  is  recognized  that  they  are  not 
independent.  At  close  distances  to  the  causative  fault. 


44 


magnitude  is  not  a  controlling  factor  for  ground  motion  in  events 
of  sufficient  size  with  fault  breakage  over  tens  of  kilometers  in 
length.  Local  effects  dominate  at  these  close  distances.  It  is 
at  greater  distances  that  magnitude  assumes  the  usual  significant 
role.  However,  the  researcher  must  exercise  care  in  treatment  of 
this  phenomena  since  an  untendad>le  position  might  occur  that  data 
could  be  extrapolated  to  show  magnitude  7.0  events  at  close 
distances  might  produce  lower  ground  motions  than  magnitude  6.0 
events  at  the  same  distance. 

It  is  significant  that  each  of  the  researchers  cited 
considers  the  uncertainty  in  the  data  and  allows  for  computation 
of  the  standard  deviation  of  the  acceleration.  This  must  be 
incorporated  in  risk  analysis  procedure  and  has  been  done  in  the 
computer  programs. 

A  topic  of  discussion  is  the  significance  of  peaks  in 
accelerograms.  It  may  be  argued  that  these  have  little  effect  on 
the  responses  of  a  structure  and  that  some  measure  of  an 
effective  acceleration  should  be  used.  Some  suggest  use  of  0.65 
times  the  peak  represents  an  effective  level.  Others  choose  a 
more  complex  approach  taking  levels  of  motion  which  compose  90% 
levels  of  motion  based  on  histogram  distribution  of  all  peaks  in 
a  record.  Although  such  approaches  are  indeed  based  on  fact  that 
structures  do  not  respond  significantly  to  single  peaks,  they 
tend  to  ignore  the  usual  engineering  approach  which  utilizes  a 
ground  motion  level  to  define  a  spectra.  The  spectra  performs 
the  function  of  giving  an  effective  level  of  motion  in  the 
frequency  response  region  of  the  structure.  It  is  thus  suggested 
that  reduction  of  ground  motion  levels  for  use  in  creation  of  a 
spectra  should  not  be  done  but  rather  the  spectra  will 
automatically  perform  this  task.  The  risk  analysis  procedure  in 
the  program  makes  use  of  mean  response  data  and  the  uncertainties 
such  as  occur  in  nature.  To  reduce  the  reported  levels  would  be 
unconservative.  It  is  substantially  better  to  utilize  a  series 
of  spectra  from  earthquake  records  to  "average"  effective 
acceleration  levels  rather  than  arbitrarily  "throw  away"  a 
portion  of  the  record.  This  is  particularly  significant  since 
the  engineering  practice  for  scaling  spectral  values  is  in  terms 
of  peak  values  and  until  some  other  value  is  fabricated  and 
agreed  upon  it  will  remain  the  only  value  clearly  established  to 
define  the  record.  The  engineer  tasked  with  design  applications 
must  exercise  caution  in  listening  to  the  academics  postulating 
improved  techniques  and  not  offering  a  corresponding  basis  for 
utilization  of  the  existing  spectral  data  base.  Without  a  data 
base  converted  to  some  measured  effective  acceleration  level  the 
engineer  would  not  be  able  to  scale  spectra  and  would  be  severely 
limited  in  characterization  of  a  site. 

Data  reported  by  Chung  (1978)  was  used  in  a  study  to  compare 
computed  and  measured  accelerations  with  distance.  Records  on 
buildings  and  vertical  records  were  excluded.  As  can  be  seen, 
there  is  significant  scatter.  The  standard  deviation  is: 
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Equation 


Mean 

Deviation* 


Standard 

Deviation 


McGuire 

Trifunac  and  Brady 
Campbell  (West) 
Ceunpbell  (East) 
Donovan  and  Bomstein 
Joyner  and  Boore 


-.005 

.089 

+  .008 

.105 

-.059 

.111 

-.041 

.101 

-.042 

.104 

-.056 

.111 

*  Minus  «  underestimates  data. 


The  McGuire  equation  appears  to  have  the  best  fit  and  least 
bias  of  the  data.  The  analysis  of  the  variances  shows  that  there 
is  a  90%  confidence  that  the  McGuire  relationship  shows  a 
significant  difference  from  the  rest,  that  is  that  within  the 
accuracy  of  the  data  the  differences  in  equations  is  meaningful. 
However,  strictly  speaking,  this  comparison  is  not  exact  because 
differences  exist  in  the  definition  of  distances  and  magnitudes. 
Chang  uses  epicentral  distance.  This  is  not  si^ificant  at  large 
distances.  In  an  attempt  to  determine  the  significance  and 
sensitivity  to  definitions  of  distance  and  magnitude,  the  Joyner- 
Boore  data  base  of  183  records  was  used  for  comparison  (based  on 
closest  distance) . 


The  data  set  used  by  Joyner  and  Boore  was  also  studied.  The 
mean  and  standard  deviation  of  the  equations  are: 


Equation 

McGuire 

Trifunac  and  Brady 
Campbell  (West) 

Campbell  (East) 

Donovan  and  Bomstein 
Joyner  and  Boore 

*  Minus  »  underestimates  data. 


Mean 

Standard 

Deviation* 

Deviatioi 

+  .001 

.0901 

+  .076 

.1885 

-.027 

.0936 

-.013 

.0913 

-.014 

.0844 

-.021 

.0830 

The  McGuire  equation  has  the  least  bias.  The  standard 
deviations  except  for  the  Trifunac  and  Brady  equations  are 
statistically  about  the  same.  It  is  not  surprising  that  the 
Joyner  and  Boore  equation  fits  this  data  best  since  it  was 
formulated  from  the  data. 


The  significance  of  selection  of  the  earthqpiake  attenuation 
equation  is  primarily  in  the  determination  of  separation 
distances.  McGuire  and  Donovan  for  example  define  the 
relationship  in  hypocentral  distance,  Trifunac  and  Brady  in  terms 
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of  the  epicentral  distance  and  Campbell  and  Joyner  and  Boore  in 
terms  of  the  closest  distance  to  the  fault.  For  risk  analysis  of 
faults  where  the  site  tends  to  be  located  near  one  end  of  the 
fault  the  selection  of  a  relationship  based  on  the  shortest 
distance  will  tend  to  give  higher  values.  Caution  should  be  used 
in  selecting  equations  like  Joyner  and  Boore  or  Campbell  that 
utilize  the  shortest  distance  for  use  in  regional  studies 
involving  floating  earthquake  epicenters  (and  not  least  distance 
to  site) .  These  relationships  will  report  significantly  lower 
accelerations  that  are  not  correct.  To  illustrate  the 
significance  of  the  attenuation  relationship  a  study  was  made  of 
a  site  0.5-miles  from  the  Hayward  fault.  Results  for  the  225- 
year  return  time  are: 


Joyner  and  Boore 

.535g 

Campbell 

.540g 

McGuire 

.396g 

Donovan  and  Bomstein 

.382g 

The  first  two  use  shortest  distance;  the  last  2  use 
hypocentral  distance.  The  hypocentral  distance  is  thought  to  be 
a  better  representation.  This  example  presents  extreme 
conditions  not  usually  encountered. 

A  study  of  the  Monterey  area  was  made.  Six  faults  were 
considered;  results  are  given  as  follows  for  the  225-year  return 
time  site  acceleration  composite  risk  from  6  faults: 

McGuire  0.2802g 

Joyner  and  Boore  0.3404g 

Campbell  0.2832g 

For  comparison  Thenhaus,  et  al.  computes  for  the  same  site 
100-year  return  time  0.20g 

500-year  return  time  0.50  to  0.62g 

The  following  projections  can  be  made.  The  100-year  return 
time  acceleration,  0.2g,  has  a  61%  probability  of  not  being 
exceeded  in  50  years.  This  translates  into  a  ratio  of  0.70  of 
the  250-year  return  time  acceleration.  The  projected  250-year 
acceleration  would  be  .02/. 07  or  0.286g.  The  500  year  return 
time  acceleration,  0.5  -  0.62g,  has  a  90%  probability  of  not 
being  exceeded  in  50  years.  This  translates  into  a  ratio  of  2.0 
of  the  250-year  return  time  acceleration.  The  projected  250-year 
acceleration  would  be  0. 5/2.0  to  0.62/2.0  or  0.25  -  0.31g.  These 
results  agree  favorably  with  those  computed  by  the  program 
(McGuire)  and  suggest  a  value  of  about  0.28g.  It  can  be  seen 
that  the  Joyner  and  Boore  value  of  0.34,  based  on  closest 
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distance  of  fault  to  site,  produces  more  conservative  results. 

To  summarize,  all  the  eUsove  cited  equations  have  been 
implemented  into  the  program  in  a  correct  consistent  manner  as 
each  author  intended  taking  account  of  each  authors  definitions. 
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PR06RAK  LIMITATIONS 


The  EARTHQUAKE  SELECTION  routine  can  process  20,000 
epicenters  to  the  file  EQUAKES. 

The  REGIONAL  STUDY  can  input  20,000  epicenters. 

The  FAULTS  STUDY  can  input  20,000  epicenters  and  30  faults. 
Each  fault  can  have  up  to  5,000  epicenters  associated  with  it. 

To  keep  the  computational  process  manageable  within  the 
limits  of  a  desktop  computer  tradeoffs  had  to  be  made  in  Monte 
Carlo  process.  As  can  be  seen  the  95  percent  confidence  bounds 
increase  as  the  probability  of  not  exceeding  the  acceleration 
increase.  This  program  is  intended  to  give  an  estimation  of 
earthquake  motions  up  to  the  90  percent  probability  of  not  being 
exceeded  in  50  years  with  high  confidence.  This  prograun  is  not 
intended  to  be  used  to  predict  5,000  year  or  10,000  year  events. 

The  user  constructs  a  model  of  the  seismicity  of  a  region  to 
which  epiceters  are  assigned,  recurrence  relationships  computed 
and  ground  motion  probabilities  determined.  It  should  be  obvious 
that  if  the  user  fails  to  include  sufficient  fault  definition  the 
model  will  lack  those  details  and  their  contribution  to  the 
site's  expected  motion.  Regional  studies  can  not  be  substituted 
for  fault  studies  with  the  same  level  of  accuracy.  While  this  is 
a  simple  but  accurate  analytical  tool  in  the  hands  of  a  trained 
engineer  it  can  be  misused  by  the  uninformed. 

WARNING 

CHECK  YOUR  RESULTS  TO  INSURE  LESS  THAN  20,000  EVENTS  WERE  SELECT¬ 
ED. 

The  program  will  display  the  number  of  events  selected  and  the 
nvunber  of  events  read.  If  5000  events  were  read  reduce  your 
area  if  possible  or  break  the  analysis  into  2  parts.  Failure  to 
do  this  will  result  in  the  omission  of  all  events  beyond  the  5000 
event  limit  thus  reducing  the  coverage. 

FOR  THE  ADVANCED  USER 

Data  and  Output  Files 

The  following  shows  the  input  data  files  created  and  the 
output  results  files  for  each  section  of  the  program. 


INPUT 


OUTPUT 


CONFIGURATION 


DEVICE . CFG 


EARTHQUAKE  SELECTION 


ONE. IN 
EPICENTR.LST 
epicenters  . EPC 


ONE. OUT 
ONE . PLT 
EQUAKES 
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REGIONAL  STUDY 


TWO. IN 
EQUAKES 


TWO. OUT 
TWO.PLT 


FAULTS  STUDY  THREE. IN  THREE. OUT 

FAULTS. LST  THREE. PLT 

Files  with  the  extension  .EPC  are  epicenter  files.  Files 
with  the  extension  .IN  are  data  files  created  by  the  program. 
Files  with  the  extension  .LST  are  files  supplied  to  the  user 
which  the  user  can  modify;  these  will  be  discussed  later.  Files 
with  the  extension  .OUT  are  output  files  to  be  printed.  Files 
with  the  extension  .PLT  are  plot  files. 

When  files  are  saved  the  above  files  are  copied  to  disk  with 
a  new  name  specified  by  the  user  as  follows: 

EARTHQUAKE  SELECTION 

ONE . IN  XXXXXXXX . IIN 

ONE. OUT  XXXXXXXX. lOU 

ONE . PLT  XXXXXXXX . IPT 

REGIONAL  STUDY 

TWO . IN  XXXXXXXX . 2 IN 

TWO. OUT  XXXXXXXX. 20U 

TWO.PLT  XXXXXXXX. 2PT 

FAULTS  STUDY 

THREE. IN  XXXXXXXX. 3IN 

THREE. OUT  XXXXXXXX. 30U 

THREE . PLT  XXXXXXXX . 3  PT 

The  program  is  set  to  use  the  standard  names  in  the  left  column. 
To  restore  a  saved  data,  output  or  plot  file  for  use  copy  the 
file  to  the  standard  name  as  follows: 

COPY  XXXXXXXX. IIN  C;\SEISMIC\ONE. IN 

COPY  XXXXXXXX. IPT  ONE. PLT 

Note  that  the  INPUT  files  are  stored  in  the  PROGRAM  DIRECTORY 
and  the  OUTPUT  and  PLOT  files  are  stored  in  the  DATA/RESULTS 
DIRECTORY.  The  program  will  re-run  or  revise  the  last  case 
executed  with  the  standard  neunes.  To  revise  or  re-run  a  case, 
only  the  INPUT  file  (.IN)  need  be  copied  to  the  standard  name. 

Epicenter  Files 

The  epicenter  data  base  is  broken  down  into  separate  files 
covering  small  regions  to  minimize  search  time.  The  user  may  add 
to  the  list  of  epicenter  file  by  editing  the  file  EPICENTR.LST 
which  contains  the  number  of  files  and  for  each  epicenter  file 
in  the  data  base  the  file  description  and  file  name. 

Number  of  Epicenter  Files 
Description  ,  xxxxxxx.EPC 
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The  epicenter  data  files  contain  data  in  the  following 
format: 


Year 

Col  5 

to 

Col 

8 

F4.0 

Latitude 

Col  20 

to 

Col 

24 

F5.3 

Longitude 

Col  26 

to 

Col 

31 

F6.3 

BM  Body  Magnitude 

Col  36 

to 

Col 

38 

FI.  2 

SM  Surface  Magnitude 

Col  54 

to 

Col 

55 

Fl.l 

OM  Other  Magnitude 

Col  61 

to 

Col 

63 

FI. 2 

IK  Local  Magnitude 

Col  83 

to 

Col 

85 

FI. 2 

Note  no  decimal  points  used. 
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Seismic  Hazard  Analysis 
Appendix 


Fault  Data  Base 


THE  FAULTS. LST  PILE 


The  FAULTS. LST  file  contains  a  list  of  over  60  faults  in 
California  whose  location  and  recurrence  have  been  defined.  This 
list  is  meant  as  a  starting  point  and  is  by  no  means  complete. 

The  values  used  are  based  on  seismicity  studies  conducted  of 
various  regions  of  California.  The  values  are  best  estimates 
based  on  the  data  available  at  the  time  they  were  made.  Program 
users  are  encouraged  to  review  the  data  and  change  the  values 
based  on  more  current  data.  To  edit  the  data  select  the  EDIT 
FAULTS. LST  option  from  the  FAULTS  STUDY  menu. 

The  user  may  advance  through  the  list  by  pressing  the  ESCAPE 
KEY;  once  the  fault  to  be  revised  is  foiind,  press  ENTER  and  move 
through  the  data  using  the  ARROW  KEYS;  enter  the  revised  numbers 
and  proceed  through  the  list.  The  program  writes  the  old  version 
of  the  list  to  a  file  labeled  FAULTS. BAK  as  a  precaution. 

The  remainder  of  this  appendix  gives  a  procedure  to 
calculate  the  Richer  recurrence  coefficients,  A  and  B  from  the 
recurrence  plots.  Note  the  program  will  attempt  to  fit  the 
recurrence  data  calculated  from  the  fault  events  as  best  as 
possible.  The  recurrence  line  should  pass  through  the  linear 
portion  of  the  data  and  usually  have  a  B  value,  slope  of  the 
line,  on  the  order  of  -.90,  the  minus  indicating  downward 
direction.  Use  this  as  a  guide  in  evaluating  the  computer 
results.  The  computer  generated  recurrence  line  is  controlled  by 
the  minimum  magnitude  allowed  in  the  analysis  and  the  distance 
parameter  which  controls  which  and  how  many  events  to  include  in 
the  calculation  of  recurrence  for  that  fault.  When  geologic  data 
is  available,  plot  the  data  on  the  recurrence  plot.  Then  move 
the  line  upward  at  the  proper  slope  to  include  the  geologic  data 
points  to  increase  fault  recurrence  if  the  geologic  data  so 
indicates . 

The  program  gives  the  user  a  great  deal  of  flexibility  in 
producing  a  valid  analysis.  Care  must  be  exercised  to  develop 
the  correct  seismic  model  of  the  earthquake  activity  of  a  region. 

Format  of  FAULTS. LST 

The  format  for  the  file  is  as  follows 

Number  of  Faults 

For  each  fault  Line  1 
Fault  Name 

For  each  fault  Line  2 

Maximum  Magnitude,  Cha’racteristic  Magnitude,  Return  time 
Hypocenter  depth  (miles) 

FAULT  TYPE  (1  ,2  or  3) 

Longitude  Point  1 


A-2 


Latitude  Point  1 
Longitude  Point  2 
Latitude  Point  2 
Longitude  Point  3 
Latitude  Point  3 

Distance  From  Fault  (miles)  for  Type  1  or  A  value 
for  Type  3 

0  for  type  1  or  B  value  for  Type  3 

Characteristic  Magnitude 

Return  Time  for  Characteristic  Magnitude 


Line  3  If  FAULT  TYPE  is  2  Define  Box  clockwise  from  upper  right 

Longitude  Point  1 
Latitude  Point  1 
Longitude  Point  2 
Latitude  Point  2 
Longitude  Point  3 
Latitude  Point  3 
Longitude  Point  4 
Latitude  Point  4 


A-3 


List  of  Faults  contained  In  FAULTS 
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Fault  -  BANNING  HaxiBUB  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.21  ,  34  PT2  116.79  ,  33.92  PT3  116.39  ,  33.83 


Figure  A-5  •  Fault  reourreuoe. 


0 


Fault  -  BIG  PINE  Maximum  Magnitude  7.6 

Fault  Longitude  /  Latitude  Coordinates 
PTl  119.85  ,  34.67  PT2  119.36  ,  34.68  PT3  119  ,  34.8 


Figure  A- 6  ,  Fault  reourrenoe 


Fault  -  CAMPROCK  Maximum  Magnitude  7 


Fault  Longitude  /  Latitude  Coordinates 
PTl  116.82  ,  34.77  PT2  116.7  ,  34.68  PT3  116.55 


34.55 


Figure  A-8  •  Fault  recurrence. 
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mx  of  Bvonta  por  Yoar  aroater  Than  or  Equal  to  Magnitude 


I  .0  1.0  2.n  3.0  U.O  5.0  6.0  7.0  6.0  9.0 

"  Magnitude 

Fault  -  CALAVERAS  Maxiaun  Magnitude  7.5 


Fault  Longitude  /  Latitude  Coordinates 
PTl  121.3  ,  36.75  PT2  121.55  ,  37.14 


PT3  122 


37.8 


Figure  A-9 


Fault  reourrenoe. 


Fault  -  CALICO  Maximum  Magnitude  7.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  116.75  ,  34.9  PT2  116.52  ,  34.65  PT3  116.22  ,  34.37 

Figure  A-io  ,  Fault  reourronoe. 


A-15 


Fault  -  COYOTE  CREEK  Maximum  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
?T1  116.52  ,  33.46  PT2  116.27  ,  33.26  PT3  115.99  ,  33.02 


Figure  A-ii  .  Fault  reourrenoe. 
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of  Bvonts  por  Yoar  Qraator  Than  or  Equal  to  Magnitude 


.0  1.0  2.0  3.0  U.O  5.0  6.0  7.0  8.0  9.0 

Magnitude 


Fault  -  DEATH  VALLEY  Maximum  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117  ,  36.53  PT2  116.84  ,  36  PT3  116.3  ,  35.6 


Figure  A-12  .  Fault  reourrenoe. 


A-17 


Fault  -  ELSINORE  Maximum  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.71  ,  33.97  PT2  116.82  ,  33.25  PT3  116.02  ,  32.8 


Figure  a-13  .  Fault  reourreaoe 


Fault  -  EMERSON  Maximum  Magnitude  7.2 

Fault  Longitude  /  Latitude  Coordinates 
PTl  116.54  ,  34.55  PT2  116.4  ,  34.44  PT3  116.26  ,  34.28 


Figure  A«i4  ,  Fault  reourronee. 
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Fault  -  fwhacE  creek  Haximu.  Magnitude  8.25 

J!?“^5,J"?3it:ude  /  latitude  Coordinates 

PTl  118.1  ,  37.64  PTl  116.85  ,  36.51  PT3  116.73  ,  35.98 


Figure  a-15  ,  reult  reourrenee. 


A-20 


Fault  -  GARLOCK  MaxiauB  Magnitude  7.75 

Fault  Longitude  /  Latitude  Coordinates 
PTl  118.9  ,  34.73  PT2  117.85  ,  35.4  PT3  116.3  ,  35.68 

rigurs  A-i« ,  Fault  raourraaoe. 


A-21 


of  Bv«nt«  p«r  y««r  Oraatar  Than  or  Equal  to  Magnituda 


.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  8.0  9.0 

Magnituda 


Fault  -  HARPER  Maxinun  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.59  ,  35.31  PT2  117.23  ,  35.1  PT3  116.92  ,  34.9 


Figure  A-ia •  Fault  raeurranoa. 
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NuBb«r  of  Bvoato  por  Yoar  Oroator  Than  or  Equal  to  Magnitude 


Fault  -  HELENOALE  Maximun  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.4  ,  34.85  PT2  117.06  ,  34.55  PT3  116.73  ,  34.25 


Figure  a*>21  .  Fault  reourrenoe. 
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0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  6.0 


S 

Fault  -  HOSGRI 

Fault  Longitude 
PTl  121.5  ,  35 


Hagnitude 


Haxinuai  Magnitude  7 

/  Latitude  Coordinates 
9  PT2  121.1  ,  35.5  PT3  120.75 


34. 


rigura  A'-22  .  Fault  raourreaoa. 


Fault  -  IMPERIAL  Maximum  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  115.57  ,  32.95  PT2  115.4  ,  32.72  PT3  115.32  ,  32.62 


Figure  a-23 .  Fault  roourrenoe 
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Fault  -  INYO  MOUNTAIN  Maximua  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.95  ,  36.92  PT2  117.81  ,  36.66  PT3  117.11  ,  36.2 


Figure  A-24 •  Fault  raourronoa. 
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Fault  -  LENWOOD  HaxinvuB  Magnitude  7.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.1  ,  34.88  PT2  116.85  ,  34.65  PT3  116.5  ,  34.3 


Vigura  A-as.  Vault  raourranoa. 
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of  Bvanto  par  Yaar  Qraatar  Than  or  Equal  to  Magnitude 


Fault  -  LIKELY  FAULT  Maxiinim  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 

121  ,  41.56  PT2  120.64  ,  41.25  PT3  120.27  ,  40.9 

Figure  a-26 •  Fault  raourranoa. 


A-31 


Fault  -  LOCKHART  Maximum  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.76  ,  35.21  PT2  117.39  ,  35.02  PT3  117.1  .  34. R« 


Figure  A-27 .  Feult  reourreaoe 


Fault  -  LUDIOW  Maximum  Magnitude  7.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  116.23  ,  34.95  PT2  116.1  ,  34.65  PT3  115.85  ,  34.45 


Figure  a<-28  •  Fault  reourreuoo. 


A-33 


Vuab«r  of  Bvonts  por  Yoar  Qroator  Than  or  Equal  to  Magnituda 


Fault  -  MENDOCINO  MaxlBvuB  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinate* 

PTl  124.45  ,  40.35  PT2  124.9  ,  40.35  PT3  125.32  ,  40.4 

Figure  a-3o •  Fault  reourreuoe. 
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Iluabar  of  Bvonta  por  Yomr  Qroator  Than  or  Equal  to  Magnitude 


Fault  -  OAKRIDGE  Naxinun  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  119.24  ,  34.23  PT2  119.05  ,  34.35  PT3  118.85  ,  34.38 


Vignre  A-33.  rault  recrarrenee. 
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Fault  -  ORTIGALITO  Maximun  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  121.34  ,  37.43  PT2  121.1  ,  37.03  PT3  120.84  ,  36.67 


Figure  a«34  •  Fault  raourreaoe 


Fault  -  OWENS  VALLEY  Maximum  Magnitude  8,25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  118  ,  36.23  PT2  118.25  ,  37.06  PT3  118.36  ,  37.6 


Flgurm  A-35 •  Fault  rmourranom. 


A-40 


Fault  '  OZENA  MaxinuiB  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  119.32  ,  34.67  PT2  119.55  ,  34.8  PT3  120.02  ,  35.04 


Flgurs  A-36 •  Fault  raourraaoa. 
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^  Magnituda 


Fault  -  PALOS  VERDE  Maximum  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  118.48  ,  33.9  PT2  118.22  ,  33.7  PT3  117.97  ,  33.33 


Figure  j^.37  .  Fault  raourrenoa. 
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Fault  -  PASAMINT  ZONE  Maxlaum  Magnitude  7.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.4  ,  36.33  PT2  117.3  ,  36  PT3  116.96  ,  35.6 

Figure  A-3«,  Fault  reourrenoe. 


A-43 


Fault  -  PINE  MOUNTAIN  Maxiaun  Magnitude  6.75 

Fault  Longitude  /  Latitude  Coordinates 
PTl  119.37  ,  34.65  PT2  119.04  ,  34.57  PT3  118.75  ,  34.62 


Figure  A-39 •  Fault  raourranoa 
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Fault  -  PIS6AH 

Fault  Longitude  /  Latitude 
PTl  116.49  ,  34.85  PT2 


Figure  a-4i 


Magnitude 


Maxlnun  Magnitude  7.25 
Coordinates 

116.23  ,  34.47  PT3  115.9  ,  34.1 

Fault  reourrenoe. 
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Iluabar  of  Bvonto  por  Tear  Greater  Than  or  Equal  to  Magnitude 


Humbar  of  Events  par  Year  Graatar  Than  or  Equal  to  Magnitude 


Fault  Longitude 


Fault  -  SAN  ANDREAS  (CRT)  Maxinun  Magnitude  8.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  118  ,  34.5  PT2  119.4  ,  34.92  PT3  121.14  ,  36.55 


Figure  A»46  •  Fault  reourreuoe. 


Fault  -  SAN  ANDREAS  (North  &  Offshore)  Maximum  Magnitude  8.25 

Fault  Longitude  /  Latitude  Coordinates 
PTl  123.65  ,  39  PT2  124.05  ,  40.03  PT3  124.9  ,  40.35 

Figure  a-48  .  Fault  reourreaoe. 
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Fault  -  SAN  CAYETANO  + 


Maxinun  Magnitude  6 . 75 


Fault  Longitude  /  Latitude  Coordinates 
PTl  119.9  ,  34.42  PT2  119.17  ,  34.45  PT3  118.85  ,  34.38 


Figure  A-49  •  Fault  reourreuoe. 


« 


Fault  -  SAN  CLEMENTE  Naxinua  Magnitude  7.7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.8  ,  32.47  PT2  118.3  ,  32.8  PT3  118.7  ,  33.19 


Figure  a-so  .  Fault  reourreaoe. 
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Nuabar  of  Events  per  Year  Greater  Tbaa  or  Equal  to  Magnitude 


Fault  -  SAN  GREGORIO  Maxinum  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  122.75  ,  37.9  PT2  122.25  ,  37  PT3  121.75  ,  36.2 

Figure  A-S2  .  Fault  reourrenoe. 
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Fault  -  SAN  JACINTO  Maximum  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.62  ,  34.32  PT2  116.9  ,  33.75  PT3  116.13  ,  33.26 


Figure  ^*53  ,  Fault  reourrenoe. 
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«  100.000 
4i 


Fault  SANTA  CRUZ  ISLAND  MaxinuB  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  120.06  ,  34.09  PT2  119.61  ,  33.97  PT3  119.4  ,  33.98 


Figure  A-ss  ,  Fault  reourreuoe 


Fault  -  SANTA  YNEZ  Maxiaum  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  120.3  ,  34.54  PT2  119.42  ,  34.49  PT3  118.9  ,  34.57 


Figure  a-57  ,  Fault  raeurraBoa. 


A-62 


Fault  -  SIERRA  NEVADA  Maxiau.  Magnitude  7.5 


/  Latitude  Coordinatee 

PTl  118.05  ,  35.27  PT2  118  ,  36.15  PT3  118.48  ,  37. 


34 


Figure  A-58  .  Fault  reourrenee. 


A-63 


Fault  -  SURPRISE  MaxiBUB  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  120.15  ,  41.9  PT2  120.2  ,  41.53  PT3  120  ,  41.15 


Figure  R-ci .  Fault  reourrcBoe. 
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of  Event a  par  Year  Greater  Than  or  Equal  to  Magnitude 
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Magnituds 


Fault  -  UN-NAMED  1  Maxinua  Magnitude  6.6 

Fault  Longitude  /  Latitude  Coordinates 
PTl  116.8  ,  35.33  PT2  116.62  ,  35.31  PT3  116.35  ,  35.3 


Figure  a-63  •  Fault  raourrauoa. 
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Nuabar  of  Bveata  per  Year  Greater  Than  or  Equal  to  Magnitude 


lIUBbar  of  Events  per  fear  Greater  Than  or  Equal  to  Magnitude 


.0  1.0  2.0  3.0  4.0  5.0  6.0  7.0  6.0  9.0 

Magnitude 


Fault  -  UN-NAMED  3  Maxiaua  Magnitude  7 

Fault  Longitude  /  Latitude  Coordinates 
PTl  117.27  ,  34.31  PT2  117.14  ,  34.44  PT3  116.53  ,  34. 


Figure  A-ss  .  Fault  reourranoe. 
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Fault  -  WHITTIER  Maximum  Magnitude  7.5 

Fault  Longitude  /  Latitude  Coordinates 
PTl  118.03  ,  33.98  PT2  117.7  ,  33.88  PT3  117.5  ,  33.75 

Figure  A-67  ,  Feult  reourrenoe. 
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DISTRIBOTION  LIST 


NAVFACENGCOM  /  CO,  ALEXANDRIA,  VA 

NAVFACENGCOM  CHESDIV  /  CO,  WASHINGTON,  DC 

NAVFACENGCOM  LANTDIV  /  CO,  NORFOLK,  VA 

NAVFACENGCOM  NORTHDIV  /  CO,  PHILADELPHIA,  PA 

NAVFACENGCOM  PACDIV  /  CO,  PEARL  HARBOR,  HI 

NAVFACENGCOM  SOUTHDIV  /  CO,  CHARLESTON,  SC 

NAVFACENGCOM  SOUTHWESTDIV  /  CO,  SAN  DIEGO,  CA 

NAVFACENGCOM  WESTDIV  /  CO,  SAN  BRUNO,  CA 

PWC  /  PENSACOLA,  FL 

PWC  /  GREAT  LAKES,  IL 

PWC  /  CO,  PEARL  HARBOR,  HI 

PWC  /  CO,  FPO  AP 

PWC  /  CO,  SAN  DIEGO,  CA 

PWC  /  CO,  OAKLAND,  CA 

PWC  /  CO,  FPO  AP 

PWC  /  PHILIPPINES,  FPO  AP 

PWC  /  YOKOSUKA,  FPO  AP 


